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ABSTRACT 

The diversity of Type la supernova (SN la) photometry is explored using a grid of 130 
one-dimensional models. It is shown that the observable properties of SNe la resulting from 
Chandrasekhar-mass explosions are chiefly determined by their final composition and some mea- 
sure of "mixing" in the explosion. A grid of final compositions is explored including essentially all 
combinations of ^^Ni, stable "iron" , and intermediate mass elements that result in an unbound white 
dwarf. Light curves (and in some cases spectra) are calculated for each model using two different 
approaches to the radiation transport problem. Within the resulting templates are models that pro- 
vide good photometric matches to essentially the entire range of observed SNe la. On the whole, the 
grid of models spans a wide range in i?-band peak magnitudes and decline rates, and does not obey 
a Phillips relation. In particular, models with the same mass of ^^Ni show large variations in their 
light curve decline rates. We identify and quantify the additional physical parameters responsible for 
this dispersion, and consider physically motivated "cuts" of the models that agree better with the 
Phillips relation, discussing why nature may have preferred these solutions. For example, models that 
produce a constant total mass of burned material of 1.1 ±0.1 Mq do give a crude Phillips relation, 
albeit with much scatter. If one further restricts that set to models that make 0.1 to 0.3 M© of stable 
iron and nickel isotopes, and then mix the ejecta strongly between the center and 0.8 Mq, reasonable 
agreement with the Phillips relation results, though still with considerable spread. We conclude that 
the supernovae that occur most frequently in nature are highly constrained by the Phillips relation 
and that a large part of the currently observed scatter in the relation is likely a consequence of the 
intrinsic diversity of these objects. 
Subject headings: supernovae: light curves, cosmology 



1. INTRODUCTION 

The application of Type la supernovae (SNe la) to cos- 
mological distance determination has yielded revolution- 
ary insig hts into the structure and composition of the 
universe l|Perlmutter et al.lll999t T Riess et al.l ll99aD. The 
utility of such explosions is based upon two empirical 
observations. First, that Type la supernovae are, for 
the most part, approximate standard candles, even with- 
out any corrections being applied. This probably reflects 
their common origin in the ex plosion of a white dwarf 
of sta ndard mass (see review bv lHillebrandt fc Niemeveil 
Il999f) . Second, a large part of the residual diversity in 
peak brightness can be corrected for by use of either light 
curve template fitting or observed correlations between 
light curve decline r ate and peak b rightness - the so called 
"Phillips Relation" l)PhilliDslll993t) . or "width- luminosity 
relation" (henceforth WLR)^. 

While the empirical relation between brightness and 
width has worked well for most purposes, as we move 
into an era of "precision cosmology" , one must feel in- 
creasing unease at the lack of an agreed upon standard 
model for how Type la supernovae explode and the pos- 
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sibility that evolutionary or environmental factors may 
erode accuracy at large distances. If the mass of ^^Ni 
made in a supernova is the dominant physical parame- 
ter affecting both brightness and light curve shape, what 
are the magnitude and direction of other possible pa- 
rameters of the explosion such as kinetic energy, inter- 
mediate mass element production, stable iron production 
and mixing? And if, as we shall find, not all physically 
plausible models obey a WLR, why has nature chosen 
to realize frequently a particular subset? That is, what 
does the WLR tell us, not just about cosmology, but 
about supernova models? If one can make progress on 
both these questions, then it may be possible to derive 
more rigorous and quantitative tools for using Type la 
supernovae for distance determination. That is the long 
term goal of our investigation. 

In this paper, we introduce a simple parameterized ap- 
proach to computing SN la explosions and use it to gen- 
erate a large grid of one-dimensional models. For each 
model, we calculate synthetic broadband light curves 
and examine the relationship between peak B magni- 
tude and i?-band decline rate. Such an approach allows 
us to study the physical parameters affecting the light 
curves of SNe la and the WLR. Several previous theoret- 
ical studies have performed similar investigations for a 
given set of models ( Khok hlov et aL 1993; Hoflich et aj 
1996, 1998; Pinto fc Eastmanll200iriHofiich et al.ll20da 
iMazzali et al..i2001l|) . These studies, however, were typi- 
cally confined to a certain class (or classes) of theoretical 
explosion paradigms. In this study, rather than adopt 
a specific theoretical framework, we take a general and 
exhaustive approach by simulating the full range of final 
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ejecta structures conceivably arising from the disruption 
of a Chandrasekhar-mass carbon-oxygen white dwarf. 

2. A SIMPLE MODEL FOR THE EXPLOSION 

2.1. Construction and Parameters of the Model 

Here a simple ansatz for the explosion model is in- 
troduced that turns out to work surprisingly well. It is 
assumed that all SNe la start from a common point, a 
1.38 Mq carbon oxygen white dwarf with a central den- 
sity 2 — 3 X 10^ g cm~^. Burning, which may be quite 
complicated, turns some part of that fuel into ash and 
deposits an internal energy in the dwarf equal to the dif- 
ference in nuclear binding energy between the initial and 
final compositions. This composition is imprinted on the 
white dwarf and, in our calculations, the corresponding 
change in nuclear binding energy is deposited uniformly 
throughout its mass. The white dwarf expands to infinity 
with a velocity distribution determined by that energy. 
The ^^Ni and ^^Co decay, and the time-dependent radi- 
ation transport is calculated. 

The major assumption that facilitates the calculation 
is that the energy deposition is shared globally. This is 
clearly true in spherically symmetric deflagration mod- 
els, since sound waves move much faster than the flame 
and share the overpressure created by the burning with 
the rest of the star. It is also true in detonations except 
for a thin layer near the surface. Throughout most of the 
star, the passage of a detonation wave changes the com- 
position and internal energy. Most acceleration occurs 
afterwards. After expansion of a factor of one million 
(before it can be seen), and the development of a veloc- 
ity field that must increase with Lagrangian mass, the 
results are the same as achieved by a simple composition 
swap and artificial energy deposition. Perhaps the best 
validation of the model is that it works. Calculations 
given in § 14.21 and § 14.31 show that the multi-band pho- 
tometry of models calculated this way agrees, both with 
previous detailed numerical simulations of the explosion, 
like the well-studied Model W7, and with a diverse set 
of observed supernovae. 

However, the simple model has one basic shortcoming. 
It will turn out that "mixing" - just how a given composi- 
tion is distributed in velocity - is an important parameter 
of the problem. It is this information that must eventu- 
ally be provided by a full, first-principle's calculation. 
Buried within this mixing parameter is also information 
about the possible asymmetry of the explosion. Such 
simulations must ultimately be three dimensional. How- 
ever, the present models are as "physical" as most other 
one-dimensional approximations that make assumptions 
about flame speeds, transitions to detonation, metallic- 
ity, etc. 

The dominant products of burning are assumed to be: 
i) ^^Ni; ii) "stable iron" , that is all the other nuclei in the 
iron group, chiefly ^^Fe and ^*Ni; and iii) intermediate 
mass elements - Si, S, Ar, Ca, henceforth referred to as 
"IME" . In the latter group Si and S are most abundant 
and, for making energy, most important, but Ca is im- 
portant for the spectrum. These are the parameters of 
the solution. Here IME ratios are adopted from Model 
DD4 of Wooslcv & Weaver (1994): by mass Si (53%), S 
(32%), Ar (6.2 %) and Ca (8.3%). Solar ratios llLndderFJ 
l2003n') would not have been much different: Si (58%), S 
(29%), Ar (7.7%), and Ca (5.2%). There is some physical 



motivation for this as well. At the temperatures where 
carbon and oxygen burn to silicon, 3-5 xlO^ K, quasi- 
eguilibrium favors an approximately solar abundance set 
(^ooslev ct al. 1973). Stable iron is taken in the calcu- 
lation to be '^''Fe. A further refinement in which stable 
iron is split into ^''Fe and ^*Ni could easily be done, but 
would add an additional parameter without greatly af- 
fecting the light curve. 

The relative amounts of these three sets of burning 
products reflect specific physical processes in the star 
and play a unique role in making the light curve and 
spectrum. "Stable iron" is a measure of burning at tem- 
peratures and densities so high (T z5 x 10^ K; pzlO^ 
g cm~'^) that nuclear statistical equilibrium is attained 
and accompanied by electron capture. To some extent, 
stable iron is also a function of the initial metallicity of 
the star (Ti mmes et al.l2003|) . Stable iron-group isotopes 
contribute to the opacity and explosion energy, but not 
to the later energy generation that makes the light curve. 
A higher ignition density, which might reflect a lower ac- 
cretion rate, increases the production of stable iron. The 
most natural location for stable iron is in the center of 
the supernova where the density was the highest, but iron 
will also exist at a level of about 5% (Z/Z©) by mass in 
the ^^Ni layer. This is a consequence of the ^^Ne present 
in the initial composition from helium burning, and the 
number is sensitive (linearly) to the metallicity. 

The double magic nucleus ^^Ni is always the domi- 
nant product of nucleosynthesis starting from a fuel with 
equal numbers of neutrons and protons (like ^^C and 
^®0) when nuclear statistical equilibrium is attained. Its 
abundance therefore reflects the extent of nuclear burn- 
ing above 5 x 10^ K at densities low enough that electron 
capture is negligible. Without doubt, it is the key player 
in the SN la light curve. Along with its decay product 
^®Co, it powers the light curve and also contributes to 
the opacity and explosion energy. The most natural lo- 
cation for the ^^Ni is also deep inside the star just outside 
the stable iron. However, turbulence, instabilities in the 
explosion, and perhaps delayed detonation may lead to 
some of the ^^Ni (and to some extent, stable iron) being 
"mixed" far out, perhaps even to the surface. 

IME are made when the density in the burning region 
declines to about 10^ g cm^"^. There, the heat capacity 
of the radiation field keeps the burning from going all the 
way to equilibrium. While these elements cause promi- 
nent line features in the SN la spectrum, their opacities 
and emissivities are not as important to the light curve 
as are those of nickel, cobalt and iron. Their production 
does, however, contribute to the explosion energy. The 
most natural location of IME is in the outer layers of the 
supernova, although instabilities and delayed detonation 
can lead to mixing among the layers in velocity space. 

Given an initial white dwarf, the free parameters of the 
calculation are thus Mni, Mpc, and Mime- To this list 
one must also add "mixing" , which is harder to quantify. 
The carbon-to-oxygen ratio and the white dwarf bind- 
ing energy (or equivalently ignition density) also affect 
the overall energetics of the explosion. Burning carbon 
releases a little more energy than burning oxygen and 
a more tightly bound white dwarf requires more energy 
to give a certain expansion speed. Here we assume a 
standard, near Chandrasekhar-mass white dwarf (1.38 
Mq) composed of equal amounts of carbon and oxygen 
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(50% by mass fraction of each) with an ignition density of 
2.9 X 10^ g cni~^. This corresponds to a net binding en- 
ergy (internal energy plus gravitational binding energy) 
of 4.62 X 10^" erg (Coulomb corrections are neglected in 
this study) . The ignition density cannot be much greater 
than this or the production of rare neutron-rich species 
will exceed the Gal actic inventory as represented by the 
sun llWooslevlll997r . The runaway will not ignite at den- 
sities below about 2.5 x 10^ g cm"-^ for reasonable choices 
of the electron screening function. 

The energy of the explosion is determined by the en- 
ergy released in the nuclear burning. The asymptotic 
kinetic energy of the SN la explosion is directly related 
to the amount of material burned in the explosion, in 
particular, for a starting composition of 50% C and 50% 
O, 

Ek = 1.56MNi + 1.74A/FO + 1.24AfiME - Eg + E,nt (1) 

where Eg ^ 3.35 B (here 1 B = 1 Bethc = lO^i erg) 
is the gravitational binding energy and Eint — 2.89 B, 
the internal energy of the progenitor white dwarf. If the 
composition were 30% C and 70% O, the yields would 
be about 7% less for Ni and Fe and 9% less for Si - Ca. 
A similar change in the opposite direction occurs for an 
oxygen-rich initial composition - 70% O and 30% C. 

This energy is deposited uniformly (as a certain num- 
ber of erg g~^) throughout the entire mass of the initial 
white dwarf and the composition of the layers changed 
appropriately (e.g., Fig.^for an explosion that made 0.7 
Mq of ^^Ni, 0.1 Mq of s^Fe and 0.3 Mg of Si-Ca), The 
ensuing evolution is fo llowed using the KEPLER implicit 
hydrodynamics code iWeaver. Zimmerman, fc Woosfevl 
119781 iWooslev et 311120021) . After less than a minute 
(star time), the explosion is homologously coasting and 
can be mapped into either of the radiation codes dis- 
cussed in i) l3.1l or 13.21 The hydrodynamical calculation 
also takes less than one minute on a desktop processor, 
thus allowing a large number of models with carefully 
controlled properties to be generated. 

2.2. Mixing 

"Mixing" is a general term referring to the fact that the 
ejected composition is not stratified into shells of pure 
ashes from a given burning process - here ^^Fe, ^^Ni, Si- 
Ca, and CO - but has somehow become blended in veloc- 
ity space. Mixing occurs at some level simply because the 
burning temperature is continuous. Carbon and oxygen 
don't burn only to Si-Ca or ^^Ni, but, at some tempera- 
tures, to a mixture of both. Stable iron may be present 
in the ^^Ni layer if the star had an appreciable (espe- 
cially super-solar) metallicity. The most difficult mixing 
to quantify however, is a consequence of the inherently 
multi-dimensional nature of the burning which, from 
start to finish, is Ray leigh- Taylor unstable and turbulent. 
Particular models like gravitationally confined detona- 
tion |Plcwa ct al. 2004; RoDkc. Wooslev. & Hillcbrandt 
I2OO6I) even have most of the burning starting at the sur- 
face and moving inwards. 

A great variety of mixing is, in principle, possible and 
to the extent that our results turn out to be sensitive 
to mixing, the simple ansatz employed for the explosion 
is questionable. However, there are some general rules 
that mixing should obey. First, mass and energy are 
conserved, so the dynamics of the explosion, including 
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Fig. 1. — Various parameterized mixing prescriptions applied to 
the model with 0.7 M© of 56Ni, 0.1 Mq of stable iron, and 0.3 
Mq of IME. The remainder is carbon and oxygen in equal (50%) 
amounts by mass. The lines for C and O lie on top of each other 
in the figure. Model C070101 is very mildly mixed corresponding 
to moving boxcar average with a window of 0.02 Mq applied three 
times. In Model M070103, the same star and composition has been 
mixed with the same 0.02 Mq interval 50 times. Model MC070103 
is extensively mixed inside 0.8 Mq and mildly mixed outside (see 
text). 



the density profile in the coasting configuration is prob- 
ably not affected very much. Second, the angle- averaged 
atomic weight of the ejecta probably decreases from the 
middle outwards. That is ^^Ni and stable iron are con- 
centrated more towards the center, carbon and oxygen 
on the outside, and Si-Ca, in between. Third, the neu- 
tron excess, 77 = T,{N — Z)X/A, also decreases from the 
center to the surface. 

A one-dimensional treatment of mixing may not be 
so bad, as far as the radiative transfer is concerned. 
The observed light curve is the emission integrated over 
the entire star. Angle-dependent mixing can certainly 
have consequences for individual spectral lines, but has 
a smaller effect on the photometry. This motivates the 
treatment of mixing as a parameterized one-dimensional 
process that stirs, but does not homogenize the compo- 
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sition. Here, that is accomplished by a running boxcar 
average of abundances in a certain mass interval moved 
from the center to surface of the star. In a typical mixing 
operation, starting with zone 1, the composition in the 
next 0.02 Mq outwards is homogenized. One then moves 
to zone 2 and does the same operation, and so on until 
the stellar surface is reached. If more mixing is desired, 
then either the mass interval is widened or the operation 
repeated, from center to surface, multiple times. The lat- 
ter approach was adopted here. Mildly and moderately 
mixed versions of a sample model are shown in Fig. ^ 

2.3. Model Nomenclature 

Models are named here according to their composition 
and the mixing prescription employed. Model Cxxyyzz is 
a "mildly mixed" model with xx/10 solar masses of ^^Ni, 
yy/10 solar masses of stable iron, zz/10 solar masses of 
IME (Si-Ca), and 1.38 - (xx-|-yy-|-zz)/10 solar masses of 
carbon and oxygen. Prior to mixing, unless otherwise 
mentioned, the ordering of the composition is stable iron 
(center), ^^Ni next, and (Si-Ca) farthest out, followed by 
unburned carbon and oxygen. By mild mixing (Fig. ^ 
we mean specifically that a moving interval of mass of 
0.02 M© was mixed from center to surface three times. 
"Moderately mixed" models, Mxxyyzz, which are most 
frequently employed in this paper, apply the same mov- 
ing boxcar average 50 times. "Highly mixed" models, 
MMxxyyzz, apply the same moving boxcar average 200 
times. 

A finer grid of the M-series models was also prepared 
for highlighting a restricted range of burned mass fS I7.2() . 
1.0 and 1.1 Mq; ^'^Ni masses, 0.4 to 0.9 Mg and sta- 
ble iron masses, - 0.3 Mq and a finer sampling of 
the intervals was also used. Additional models were also 
constructed to explore different mixing prescriptions de- 
scribed in § 16.2.21 and § 17.51 Figure ^ shows the compo- 
sitional structure of a representative model with several 
mixing prescriptions. 

3. RADIATION TRANSPORT CALCULATIONS 

The radiative transfer calculations, which by far con- 
sumed most of the computer time, were performed using 
two independent radiative transfer codes which employ 
very different numerical methods and atomic data. Both 
start with a supernova that has already expanded to the 
"coasting" or "homologous expansion" phase. Opera- 
tionally, that means the explosion model from kepler 
was linked into the radiation transport code at an age of 
10,000 s. 

3.1. Monte Carlo - Sedona 

The SEDONA code (|Kasen et al.l l200(Tj) is a time- 
dependent multi-dimensional Monte Carlo radiative 
transfer code, designed to calculate the light curves, 
spectra and polarization of supernova explosion mod- 
els. Given a homologously expanding SN ejecta struc- 
ture, SEDONA calculates the full time series of emergent 
spectra at high wavelength resolution. Broadband light 
curves are then constructed by convolving the synthetic 
spectrum at each time w ith the appropriate Bessel filter 
transmission functions (Bessell 1990). SEDONA includes 
a detailed treatment of gamma-ray transfer to determine 
the instantaneous energy deposition rate from radioac- 
tive ^^Ni and "^^Co decay. Radiative heating and cooling 



rates are evaluated from Monte Carlo estimators, and the 
temperature structure of the ejecta determ ined bv iter- 
ating^the model to thermal equilibrium. See lKasen et alJ 
( 2006) for a detailed code description and verification. 

Several significant approximations are made in SE- 
DONA, notably the assumption of local thermodynamic 
equilibrium (LTE) in computing the atomic level pop- 
ulations. In addition, bound-bound line transitions are 
treated using the expansion opacity formalism (imply- 
ing the Sobolev approximation). Although the SEDONA 
code is capable of a direct Monte Carlo treatment of 
NLTE line processes, due to computational constraints 
this functionality is not exploited here. Instead, the line 
source functions are treated using an approximate two- 
level atom approach (see Eq. 12 next section). In the 
present calculations, we assume for simplicity that all 
lines are "purely absorptive", i.e., in the two-level atom 
formalism the ratio of the probability of redistribution to 
pure scattering is taken to be eth — 1 for all lines. In this 
case, the line source functions are given by the Planck 
function, consistent with our adoption LTE level popu- 
lations. The one exception is the calcium lines, which are 
assumed to be p ure scattering (eth = 0) for the reasons 
discussed in iKa scn (2006). 

The assumption of purely absorbing lines has 
been common in pre vious SN tra nsfer calcula - 
tions llPinto fc Eastman 2001; Blinniko veTall |2006). 
iKasen et alJ 1)2006..') demonstrate that the eth = 1 ap- 
proach does in fact capture the true NLTE line fluores- 
cence processes operative in SNe la light curves, however 
quantitative errors in the broadband magnitudes are ex- 
pected on the order of 0.1 to 0.3 mag. A refined calibra- 
tion would better represent the true NLTE redistribution 
probabi lities and may provide more accurate results (e.g., 
lHoflichl[T995|l . In particular. fKasen et all l)2006(l shows 
that for SN la models assuming eth < 1 leads to more ac- 
curate model colors in the maximum light spectrum. For 
this reason, our model peak magnitudes, colors and de- 
cline rates should be considered uncertain by this moder- 
ate amount. The adoption of an alternative value eth 7^ 1 
would lead to a shift in the location of the models in the 
WLR and color plots discussed below, however this shift 
is most likely in a uncorrelated way that does not sig- 
nificantly change the slope of the model relation or the 
level of dispersion. On the other hand, if eth depends in 
a systematic way on temperature or density (as to some 
extent it must) this could in principle lead to correlated 
errors that affect the slope of the model relation. The 
issue will be the subject of future investigations. 

The two-level atom framework applied here is just one 
of several uncertainties that affect the radiative transfer 
calculations. In addition, inaccuracy or incompleteness 
in the atomic line data can be a source of significant er- 
ror (see S I4.1|I . The inaccuracy of the LTE ionization as- 
sumption may also h ave signi ficant consequences for the 
B-band light curves l|Kasen fc Wooslcv 2006). At later 
times (> 30 days after i?-band maximum) the NLTE ef- 
fects become increasingly significant and the model cal- 
culations become unreliable. One should keep in mind 
all the sources of uncertainty when considering the model 
relations discussed below. 

The numerical griding in the present calculations was 
as follows: spatial: 120 equally spaced radial zones with 
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Fig. 2. — Logarithm of the density profile of unmixed Model 
C070203 computed by STELLA at t = 90 days (red lines) compared 
to the structure obtained by homologous expansion from the initial 
model (black lines). 



a maximum velocity of 30,000 kms^^; temporal: 100 
time points beginning at day 2 and extending to day 80 
with logarithmic spacing A logt = 0.175; wavelength: 
covering the range 100-30000 A with resolution of 10 A. 
Extensive testing confirms the adequacy of this griding 
for the problem at hand. Atomic line list data w as taken 
from the Kurucz CD 1 line list l|Kurucz 1 119941) . which 
contains nearly 42 million lines. 

A total of 10^ photon packets were used for each 
calculation, which allowed for acceptable signal-to-noise 
(S/N) in the synthetic broadband light curves. For a 
few models, further time-independent calculations were 
performed using additional photon packets in order to 
calculate high S/N synthetic spectra at select epochs. 

3.2. Multi-energy group diffusion - Stella 

The photometry of most of the models, including all 
of them with the standard "M-mixing" , was also in- 
dependently verified using an ol der multi-energy group 
radiation hydro code Stella ijBlinnikov et al.l Il998t 
^linnikov & Sorokina 2000). Stella was developed 
primarily for Type II supernovae, where the effects 
of coupling of radiation transfer to the hydrodynam- 
ics are more imp ortan t, e.g., during shock propagation 
l|Blinnikov et alJl2003t fchugai et al1l2004() . Stella is 
not a Monte Carlo code, but employs a direct numerical 
solution of radiative transfer equation in the moment ap- 
proximation. Time-dependent equations for the angular 
moments of intensity in fixed frequency bins are coupled 
to the Lagrangian hydro-equations and solved implicitly. 
The photon energy distribution may be quite arbitrary. 
Due to coupling with hydrodynamics, which is generally 
not needed for the coasting phase of SN la (though see 
below), the radiative transfer part of the calculation is 



somewhat cruder than in SEDONA. 

The effect of line opacity is treated by Stella, in 
the current work, as an expansion opacity according to 
lEastman fc Pintol (|199^, similar to SEDONA. The line 
list is limited to ~160 thousand entries, selected from 
the strongest down to the weakest lines until saturation 
in the expression for expansion flux opacity is achieved. 
This is, in general, much less than in SEDONA, moreover, 
the list is optimized for solar abundance only. 

The ionization and atomic level populations are de- 
scribed by Saha-Boltzmann expressions. However, the 
source function is not in complete LTE. The source func- 
tion at wavelength A is 

= ethSA + (l-eth) Ja, (2) 

where eth is the thermalization parameter, B\ the Planck 
function, and J\ the angle mean intensity. In this work, 
STELLA assumes eth = 1- Hydrodynamics coupled to ra- 
diation is fully computed (homologous expansion is not 
assumed). This question, and major other recent im- 
pro vements in the code St ella (introduced after the pa- 
per Blinnikov et al. 1998', was published) are described 
in (Blin nikov et al. 200,^). 

The heating by the decays of ^^Ni ^'^Co^ ^^Fe is 
taken into account. It is assumed that positrons, born 
in the decays, are trapped so they deposit their kinetic 
energy locally. To find the radioactive energy deposition, 
we treat the gamma-ray opacity as a pure absorptive 
one, and solve the gamma-ray t ransfer equation in a one- 
group approximation following iSwartz et al.l l) 19951) . 

The effective opacity is assumed to be — 
0.05y"e cm^/g, where Ye is the total electron number 
density divided by the baryon density. Although this 
approach is checked against other algorithm s, e.g. those 
used in eddington ijBlinnikov et al.llT99^ it is less ac- 
curate than the full Monte-Carlo treatment in SEDONA 

To calculate SNe la light curves Stella can use up 
to 200 frequency bins and up to '^400 zones in mass as 
a Lagrangian coordinate on a modest processor, but all 
current results are obtained with 100 groups in energy 
and 90 radial mesh zones. 

As described by Blinniko v et all 1)19981) . Stella has 
only an approximate treatment of light travel time cor- 
rection, since it works not with time-dependent intensity, 
but with time-dependent energy and fiuxes only. SE- 
DONA is superior in this aspect since it works directly 
with packets of photons. 

The approximation of homologous expansion is usu- 
ally exploited in the radiative transfer codes that 
neglect h ydrodynamics ifEastman fc Pintol [19 93s; 'Lucd 
1200 5: Kas en et al.ll200d) However. iPinto fc Eastman! 
l)2000fl point out that the energy released in the ^^Ni 
decay can infiuence the dynamics of the expansion. The 
^^Ni decay energy is 3 x 10^^ erg g~^, which is equivalent 
to a speed of 2.5 x 10'^ km s~^, if tra nsformed into the 
kinetic energy of a gram of material. IPinto fc Eastmanl 
l)2000(l state: "Since the observed expansion velocity of 
SNe la is in excess of 10* km s~ , we expect that this 
additional source of energy will have a modest, but per- 
haps not completely negligible, effect upon the velocity 
structure." In reality, the heat released by the ^^Ni de- 
cay will not all go into the expansion of the SN la. If 
the majority of ^^Ni is located in the central regions of 
the ejecta then the main effect is an increase in the en- 
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Fig. 3. — Comparison of the UBVRI light curves calculated by 
STELLA (dashed line) and SEDONA (solid line) for Models M070103 
(left) and M040303 (right) Agreement for the _B-band light curve 
during the first 20 days after maximum, which is of greatest interest 
in this paper, is excellent for the models that produce medium to 
large quantities of ^^Ni and still quite good for those with low 
values like M040303. The divergence in the models after day 25 
is due primarily to the lesser amount of lines used by STELLA to 
compute the expansion opacity, and may also be influenced by the 
different treatment of the Ca II IR triplet lines i Kascn 2009) . 



tropy and local pressure (both quantities are dominated 
by photons in the ejecta for the first several weeks). The 
weak overpressure will lead to a small decrease in den- 
sity at the location of the 'nickel bubble', as well as to 
some acceleration of matter outside the bubble. It is of- 
ten said that the expansion of the ejecta is supersonic 
and that pressure cannot change the velocity of the mat- 
ter, but one should remember that in the vicinity of each 
material point we have a 'Hubble' flow, so differential 
velocities are in fact subsonic in a finite volume around 
each point. 

All of the models considered here have moderately 
mixed distributions of ^^Ni. Other models (without mix- 
ing) demonstrate that the nickel bubble, i.e. the depres- 
sion of density in ^^Ni-rich layers, continues to grow dur- 
ing the coasting stage. The effect is modest, but as ex- 
pected it may result in a ~10% difference in velocity and 
density This effect is evident, e.g., for Model C070203 
(see Fig. The solid black line is the initial model 
scaled to our result at 90 days since explosion. We see 
that the density profile has changed due to the ^^Ni and 
^^Co decays; it clearly deviates from homology. 

The change in the density is important for the deposi- 
tion of gamma-ray energy, which is reflected in the light 
curve. This may explain partly some of the differences 
between Stella and SEDONA results. The larger the ^^Ni 
mass, the larger is the difference expected. 

4. CODE AND MODEL VERIFICATION AND VALIDATION 

4.1. Comparison of Results from Sedona and Stella 

Figure |31 shows the comparison of the UBVRI light 
curves calculated by STELLA and SEDONA for Models 
M070103 and M040303. Agreement for the B-band hght 
curve during the first 20 days after maximum, which 
is of greatest interest in this paper, is quite good for 
the models that produce medium to large quantities of 



^®Ni, and still reasonable for those with low values like 
M040303. In general, the results of the two codes are 
mutually confirming, however one does note moderate 
differences in the _B-band light curve decline rates and 
peak magnitudes, and substantial differences in the early 
/ and i?-band light curves and the late time B and V- 
band light curves. Note that the results of both codes 
become increasingly unreliable after 40 days past B-band 
maximum, due to the increasing inadequacy of the LTE 
excitation /ionization assumption. 

Our numerical experiments suggest that these differ- 
ences in the radiative transfer results are explained pri- 
marily by the different atom ic line data u sed in the sep- 
arate codes. As discussed in iKagen! |200^ , the use of an 
extensive atomic line list (with > 5 million lines) is crit- 
ical in synthesizing the light curves of SNe la, especially 
in the red and near- infrared wavelength bands. Because 
the number of weak lines treated in SEDONA is much 
larger than that in Stella, the results of the former code 
show a generally superior correspondence with observa- 
tions. However, even the more extensive atomic line list 
used in SEDONA is likely still somewhat incomplete and 
inaccurate. Inadequacies in the available atomic data re- 
main an important source of uncertainty in supernova 
light curve modeling. 

Other differences between the two transfer codes may 
also contribute to the discrepancies seen in Figure |3| SE- 
DONA employs a finer frequency grid, and thus better re- 
solves individual spectral features. The features can have 
an important impact on the broadband magnitudes. In 
addition, SEDONA treats the Ca II IR triplet lines as pure 
scattering, which likely better represents the source func- 
tion. Because the Ca II lines become excessively strong 
at later e pochs, they c an affect the radiative transfer in 
all bands (•Kase nl2006j) . SEDONA also includes a more de- 
tailed, multi-group gamma-ray transfer procedure, and 
thus more accurately determines the radioactive energy 
deposition rate and late times. The main advantage of 
STELLA is its self-consistent treatment of hydrodynam- 
ics, but arguments in the previous section and the den- 
sity plot Fig. [3 suggests that deviations from homology 
are not large. For present purposes, SEDONA is superior 
and most of the subsequent figures and discussion in this 
paper are based on the results of that code, except where 
otherwise noted. 

4.2. Comparison With Analytic Approximations and 
Other Calculations 

The primary output of the hydrodynamic explosion 
calculation is the final density structure (e.g., density 
versus velocity) of the ejecta once it has reached the ho- 
mologous phase. Interestingly, the final density profiles 
of all models are well characterized by a simple exponen- 
tial function which depends only upon the kinetic energy 
of the explosion 



p{v,t) 



where Ve — 



6M. 



= 2455 



Mch 

S-KV^t^ 

Ek 



exp[-t;/'i;e] (3) 



ch 



1051 



1/2 



km 



(4) 



Figure 01 shows that this simple analytic formula holds 
very well for all but the innermost ejecta [v < Ve)- For 
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Fig. 4. — Logarithm of the density structure of explosion models 
with different kinetic energy (black lines) compared to the expo- 
nential function (Equation]^ red lines). The kinetic energy Efc of 
each model is marked on the figure. 
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Fig. 5. — Logarithm of the density structure at 10000 sec for 
Model M060303 (black line) compared to that of W7 (red line). 



the high Ek models, the formula overestimates the cen- 
tral densities, while from the low Ek models it underes- 
timates them. 

The properties of our models can be compared to 
existing standard SN la explosion models. A much 
studied 1-D model that has been shown to agree with 
typical observed light curves and s p ectra is Model 
W7 of iNomoto. Thielemann. & Yokoil j lg84tl (see also 
iThieleniann^tani^Sra^anK^^^ In terms 

of our four composition parameters, the abundances in 
the final frozen-out Model W7 are ^"^Ni, 0.59 M©, stable 
iron (54.56pe, 55Mn, ss^eojs^j)^ q 35 Mq, and Si-Ca, 0.27 
Mq. If our simple model for the hydrodynamics is cor- 
rect and if mixing is not a major issue, this should be a 
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Fig. 6.— Demonstration that the simphfied Model M060303 
(black lines) well reproduces the observable properties of the well- 
known W7 model (red lines). Left: Comparison of the near max- 
imum light (texp = 18 day) spectra of the two models. Right: 
Comparison of the SVR-band light curves of the two models. 



good match for our Model M060303 (which has 0.6, 0.3, 
and 0.3 M© of ^^Ni, stable iron and IME respectively). 

In Fig. El we compare the density profile of Model 
M060303 to that of W7. The agreement is quite good 
in the inner layers, although there are discrepancies for 
velocities v > 13,000 km s~^. This is the region in 
which the deflagration burning began to be quenched 
in W7, which lead to the production of a density spike. 
This density spike would probably be absent in a multi- 
dimensional simulation. 

In Fig. |H| we compare the light curve and near max- 
imum hght spectra (ioxp= 18 days) of Model M060303 
and W7. Despite the differences in the density pro- 
file and approximate representation of the composition 
in M060303, the overall agreement is reasonably good. 
The light curves of M060303 are slightly faster than W7, 
due to the slightly lowe r dens ities and greater stable iron 
group production (see ii6.2.2|l . The maximum light spec- 
tra show only minor discrepancies in the lines of calcium 
and in the ultraviolet reflecting the different mixing em- 
ployed here and the fact that the Ca abundance in Model 
M060303 is about twice as large as in W7 (0.024 Mqvs 
0.012 Mq) and extends to higher velocity. 

4.3. Comparison to Individual Observations 

Within our model grid are examples that resemble a 
broad range of observed SNe la, from photometrically 
normal events like SN 2001el, to bright, broad events 
hke SN 1991T, and faint, narrow ones like SN 1991bg. 
In Fig. El the UBVRI light curves Model M070103 
(AfNi = 0.7 M0) are compared with those of the nor- 
mal Type la SN 2001el ( Krisciunas et al. 2003). The 
observed light curves of SN 2001el have been corrected 
for dust extinction using the estimates = 0.57 and 
— 2.88, and the distance modulus is taken to be 
^I. = 31.3 ijKrisciunas et all 120031 iWang et alJ l2003j) . 
Generally good agreement between model and observa- 
tions is found in all bands. The model further reproduces 
the double-peaked behavior of the I-band light curve, al- 
though the secondary maximum is stronger than in the 
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Fig. 7. — Right: Comparison of the broadband synthetic hght 
curves of Model M080202 to observations of the bright/broad 
Type la SN 1991T (Lira et al. 1998). Left: Comparison of the 
broadband synthetic light curves of Model M010309 to observa- 
tions of the Type la S N 1999by, a subluminous SN 1991bg-like 
event IGarnavich et ar.ll2004MV 





-10 10 20 30 40 50 
Doys Since B-Bond Uoximum 



4000 5000 6000 7000 8000 9000 
Wovelength (Angstroems) 



Fig. 8. — Verification of the Model M070103 against SN la obser- 
vations. Left: Synthetic broadband light curves of the model (solid 
lines) compared to obs ervations of the normal Type la SN 2001el 
IKrisciunas et aT1l2003l filled circles). In order to improve the vi- 
sual comparison, the observations have been offset by —0.1 mag 
in order to better align them with the model light curves. Right: 
Synthetic spectra at three different epochs with respect to B-band 
maximum of the model (black lines) compared to observations of 
the normal Type la SN 1994D (red lines). 



observations. 

FigureHdemonstrates that our models also include ex- 
amples resembling more extreme SNe la events. The left 
panel of the figure compares the light curves of Model 
M080101 ( Mn\ = 0.8 Mn) to those of the overly bright 
SN 1991T l)Lira et al.lll998() . We adopt th e reddening es- 
timate of Ay = 0.50. i?„ = 3.1 from PhiUips et alJ l|1999f) 
and the Cepheid distance me asurement of ^ = 30.56 
from iGibson fc Stetsonl l)2001|) . although there is sub- 
stantial uncertainty in these values. The right panel of 
the figure compares the light curves of Model M010207 
{Mm = 0.1 Mq) to those of the subluminous SN 1991bg- 
like event SN 1999by. The reddening and distance here 
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Fig. 9. — Comparison of the broadband synthetic light curves 
of Model M040303, computed by STELLA (solid lines), to observa- 
tions of the Type la SN 1994D IRichmond et al.ll995tlMeikle etaTl 
11993) 



are taken to be A; = 0.43, Rv — 3.1, — 30.75 
(jGarnavich et alJl2004|) . Good agreement with observa- 
tions is found for both examples. 

Full time-series of synthetic spectrum calculations are 
available from the SEDONA code and Fig.|Sl(right panel) 
compares the Model M070103 spectra at several epochs 
to observations of the norma l Typ e la SNe 1994D 
l|Meikle et alJlT99l IPatat et alJlT99fil) . While there are 
some differences in detail, the model reproduces quite 
well the essential spectroscopic features and colors over 
a wide evolution time. The only major discrepancy oc- 
curs in the Ca II IR-triplet features, which is too strong 
in the model. This confirms the adequacy of the transfer 
calculations to model the spectral and color evolution of 
SNe la over the timescale of interest. 

Fig. compares synthetic light curves of Model 
M040303, compute d by Stella, to observati ons of the 
Type la SN 1994D l|Richmond et alJil99,5: .Meikle et alJ 
Il99fij) . The agreement is especially good in UBV, and 
in general it is not worse than fo r the MPA deflagration 
models presented by Blinnikov e t al.l ||200(tD . Compari- 
son with SEDONA results for the same model in Fig. |21 
suggests that the agreement with observations in R and 
/ filters can be improved by extending the line list used 
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in STELLA for computing expansion opacity. 

5. BASIC PHYSICS OF THE WIDTH-LUMINOSITY 
RELATION 

Before turning to the model calculations, it is useful to 
summarize the basic radiative transfer physics relating 
to the brightness and decline rate of SN la light curves. 
A more detailed discussion of the transfer effects can b e 
found in a companion paper ijKasen fc Wooslevll20nl . 
Because the light curves of SNe la are influenced by 
a number of physical parameters, the behavior of our 
model light curves can not usually be explained by ref- 
erencing a single cause, but rather require consideration 
of a combination of interrelated transfer effects. We at- 
tempt to describe the most important effects individually 
below. 

It is well known that the light curves of SNe la are 
powered entirely by the decay of radioactive ^^Ni (and its 
daughter ^^Co) synthesized in the explosion. The mass 
of ^^Ni produced (Mni) is therefore the primary determi- 
nate of the peak brightness of the event. On the basis of 
approximate analytic models, l)Arnett "1982'l showed that 
the bolometric luminosity at peak is roughly equal to the 
instantaneous rate of radioactive energy deposition 

Lp fMmew{-tp/tNt) (5) 

where tp is the rise time to peak, tNi ~ 8.8 days is the 
^^Ni decay time and / is the percentage of the gamma- 
ray decay energy that is trapped at the bolometric peak 
(typically /> 0.9). 

In SNe la, the ejecta remain optically thick for the first 
several months after explosion. The width of the bolo- 
metric light curve is related to the photon diffusion time. 
The basic diffusion physics can roughly be understood 
using simple scaling arguments. In a standard random 
walk, the diffusion time is given by td ^ E? /\pC where 
R is the radius, Xp = l/^cp is the photon mean free path, 
and K is the mean opacity. In homologous expansion 
R = ipwivt where v is the characteristic ejecta velocity, t 
is the time since explosion, and (jiNi is a factor describing 
the fractional distance between the bulk of ^^Ni and the 
ejecta surface, roughly: (pm ~ — M^^)/M where M^j 
is the center of mass of the ^^Ni distribution. Using the 
scaling relations for the total ejected mass M ~ pv^t^ 
and kinetic energy Ek ~ Mv'^, one finds 

td ^ 0Ni k'/'M'/''E-'/\ (6) 

The important parameters affecting the bolometric dif- 
fusion time are thus the total mass M, the kinetic 
energy Ek, the radial distribution of nickel (/)Ni, and 
the effective opacity per unit gram, n. The last of 
these is the most complicated, depending upon the com- 
position, density, and thermal state of the ejecta, as 
well as the velocity shear across the ejecta. In gen- 
eral, K increases with temperature/ionization, thus mod- 
els with larger Mni will typically have slightly longer 
diffusion times (Khokhlov et al. 1993: Pinto fc E astmanI 
1200 1: Hoflich et al. 2002; .Kasen fc Wooslev 20061 

A further factor influencing the bolometric light curve 
is the rate at which the ejecta become transparent to 
gamma-rays from radioactive decay. Near maximum 
light, the densities in a Chandrasekhar-mass model are 
high enough that nearly all gamma-rays are trapped lo- 
cally (/ > 0.9). However, by 15 days after maximum 



densities have dropped such that a substantial percent- 
age of the gamma-rays escape the ejecta without being 
thermalized. Models in which the bulk of ^^Ni is located 
further out in mass coordinates will experience a more 
rapid transition to gamma-ray transparency and hence 
possess a generally faster bolometric decline rate. 

In addition to the parameters affecting the bolometric 
decline rate just mentioned, one must also consider the 
physics affe cting the spectroscopic an d color evolution of 
SNe la. InE asen fc Wooslev! l|2006H . it was shown that 
the WLR arises primarily from a broadband effect. In 
particular, the i?-band light curve decline rate depends 
sensitively on the rate at which the SN colors shift pro- 
gressively redwards following maximum light. Dimmer 
SNe la (i.e., those with lower Mni) exhibit a generally 
faster color evolution, which is the primary reason for 
their faster i3-band decline. Physically, this reflects the 
faster ionization evolution of dimmer SNe la. Follow- 
ing maximum- light, the SN colors are increasingly deter- 
mined by the development of numerous Fe II and Co II 
lines that blanket the bluer wavelength bands and, at the 
same time, increase the emissivity at longer wavelengths. 
Because dimmer SNe la are generally cooler, they expe- 
rience an earlier onset of Fe HI to Fe II recombination in 
the iron-group rich layers of ejecta. Consequently, Fe II 
and Co II line blanketing develops more rapidly in dim- 
mer SNe la, resulting in a more rapid evolution of the 
SN spectral energy distribution to the red. This is the 
principle explanation for their faster B-band decline rate. 

As a corollary, one realizes that the velocity distribu- 
tion of iron group elements plays an additional important 
role in determining the broadband light curves. Models 
in which iron group elements are concentrated at low ve- 
locities are unable to form strong Fe 11/ Co II line features 
in the post-maximum epochs, and consequently will ex- 
hibit a slower color evolution (and hence _B-band decline 
rate) compared to those in which the iron group elements 
are mixed out to higher velocity. 

6. RESULTS 

6.1. All Models Combined 

The WLR is often quantified as relation between peak 
B-band magnitude Mb and the drop in i?-band mag- 
nitude 15 days after peak AMi^^B). Figures [TOI and 
im show the WLR for the full set of moderately mixed 
models calculated with the SEDONA and STELLA codes, 
respectively. The models are color-coded by the mass of 
^^Ni each produces as defined in Fig. EH Overplotted in 
both figures, as a sh aded band, is the empirical WLR of 
iPhillips et al.l l)1999j1 . for which an absolute calibration 
of Mb = —19.3 mag at AMx_5lj?) = 1 -1 mag and a dis- 
persion oi a ~ 0.15 mag l|Hamuv et al.lil996|1 have been 
adopted. 

In contrast to the observed behavior, the models span 
a wide region in the figure. Contrary to popular expecta- 
tions, the light curve width and luminosity are not both 
determined by a single parameter, the ^^Ni produced in 
the explosion. If this were the case, all models with a 
given mass of ^^Ni would collapse to one point on the 
plot, and that point would be in the gray band. 

As the plots confirm, A^ni is clearly the dominant pa- 
rameter affecting the SN peak magnitude. For a given 
Mni value, Mb varies only by about ±0.25 mag. The 
decline rate AA/i5(i?), on the other hand, spans at least 
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Fig. 10. — Relationship between the B-band decline rate 
AMi5(iJ) and peak _B-band magnitude for the full set of M-series 
models. The color coding shows the ^^Ni mass, which varies from 
0.1 to 1.1 M0. The shaded region is the observed width-luminosity 
relation of|Phillips et al. 1 1999) with a calibration AIb = —19.3 for 
AMi5(i?) = 1.1 and a dispersion of tr = 0.15 mag. In contrast to 
the observations, the models occupy a wide region in the plot, in- 
dicating a sensitivity of the light curves to parameters other than 
^^Ni.The systematic offset of ~ 0.06 magnitudes in color between 
the observations and models suggests that assumed thermalization 
parameter en^ = 1, likely overestimates the actual redistribution 
probability. 
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Fig. 11. — Relationship between the B-band decline rate 
AMi5(iJ) and peak B-band magnitude for the full set of M-series 
models computed by STELLA. The color coding is the same as in 
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a full magnitude for a given Mni mass, indicating its 
sensitivity to additional physical parameters. Given the 
scaling of the diffusion time (Equation EJ, one can antic- 
ipate that two very important parameters are the total 
kinetic energy of the explosion, Ek, and the radial dis- 
tribution of ^^Ni, 0Ni- This is confirmed in the following 
sections by examining suitable subsets of models. 

Careful examination of the figures shows, in fact, that 
for fixed Mni, the model light curves actually exhibit 
an "anti-Phillips relation", i.e., the brighter supernovae 
are more narrow. This is not a surprising result, as SNe 
with shorter diffusion times (i.e., faster light curves) lose 
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Fig. 12. — Relationship between the maximum light B-V color 
and the peak B-band magnitude for the full set of M-series models. 
The color coding shows the ^^Ni mass, which vari es from 0.1 to 
1.1 M (7i. The solid line is the observed relation of IPhillips et ail 
119991) . with a calibration Mb = -19.3 for AMi5(B) = 1.1. 



a smaller percentage of their internal energy to adiabatic 
expansion, and thus reach a brighter peak earlier. Equiv- 
alently, this can understood as an expression of Arnett's 
rule (EquationEl). Models with broader light curves typ- 
ically have a longer rise time tp and thus, for given Mni, 
are dimmer at peak. 

Although the models in Figs. 1101 and 1111 do not repro- 
duce the observed WLR relation, they nonetheless lead 
to a very important physical insight - the true SN la 
explosion mechanism realizes only a small subset of the 
theoretically conceivable possibilities, implying a rather 
tight internal correlation between the relevant physical 
parameters. This places a strong constraint on theoreti- 
cal explosion paradigms. Indeed, in the following we use 
this constraint to deduce some of the properties of the 
SN la ejecta structure. 

The large spread seen in Figs . 11111 and ITTl also suggests 
that intrinsic SN variation is likely a significant source of 
intrinsic scatter in the WLR. Because AA/i5(i?) depends 
upon other parameters than M-^i, any uncorrelated vari- 
ation of the secondary parameters leads to dispersion in 
the WLR. Any systematic variation of the parameters 
with progenitor environment could form a potential ba- 
sis for evolutionary effects. 

Interestingly, if one plots the peak model magnitudes 
of the entire set vs the B-V at maximum instead of 
AMi5(i3), a tighter correlation results (Fig. EJ- The 
slope of the model correlatio n closely reseinbles t hat of 
the observed relation given in IPhillips et aLlllT99l . The 
models are systematically redder than the observations 
by ^ 0.06 mag, suggesting that our assumed thermal- 
ization parameter (eth = 1) likely overestimates the true 
redistribution probability in SNe la I3.1|) . Though it 
is clear that not all the models shown in Fig. are 
frequently realized in nature, the smaller dispersion sug- 
gests that such color indicators may be less sensitive to 
intrinsic variation in the supernovae. 

6.2. Sources of Dispersion 
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Fig. 13. — Same as Fig. 1101 but color coded to show the total 
burned mass of the models. Models with the same total explosion 
energy, i.e., points of the same color do roughly yield a width- 
luminosity relation in which "brighter equals laroader" . The best 
agreement is for a burned mass of l.litO.l M0, though the narrower 
light curves are somewhat fainter than the observations if burned 
mass is a constant. 
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Fig. 14. — Effect of the total burned mass (and hence kinetic 
energy) on the B-band light curve for a subset of M-series models 
with constant ^®Ni. The models shown have MiMi= 0.5 Mq and 
AfFo= 0.1 M0 and AfiME= [0.1,0.3,0.5,0.7 Mq]. The left panels 
shows the B-band model light curves, while the right panel shows 
the corresponding width-luminosity relation (as in Fig. 1101 . For a 
constant ^'^Ni mass, models with a greater total burned mass have 
higher kinetic energy, and hence faster light curves. All else being 
equal, the full width of the observed region corresponds to a range 
in burned mass of only 0.2 Mq. 



The large spread in the model WLR of Figs.llUlandllll 

indicates that parameters other than Mni significantly 
affect the light curves of SNe la. We show below that 
the most important of these are the total burned mass, 
Mburn, which determines the Ek, and the stable iron 
mass, Afpc, which influences (pm- In addition, the de- 
gree of direct ^^Ni mixing is significant as well. For a 
given Mni these parameters have significant impact on 
the decline rate AMi5(i3), and hence may act as sources 
of dispersion in the WLR. Here the effects of each are 
quantified using suitable subsets of the models 

6.2.1. Effect of the Total Burned Mass 

The models shown in Figs.llOlandllllspan a wide range 
in the fraction of the original white dwarf which is burned 
to heavier elements. Because the nuclear energies re- 
leased in burning carbon and oxygen to iron, ^^Ni or IME 
are all about the same, and since the same initial white 
dwarf is used for all calculations, the total mass burned in 
the explosion, Mburn, essentially dictates the final kinetic 
energy, Ek, of the ejecta. Given the inverse dependence 
of the diffusion time on Ek (Equation EJ, a higher value 
of Mburn can be expected to lead to a relatively faster 
bolometric light curve for any given model. In addition, 
because the velocity of the bulk of ^®Ni ejected in the 
explosion increases with Ek, a higher Mburn should also 
lead to a faster evolution of the model colors and gamma- 
ray transparency, both of which further contribute to a 
faster B-band light curve decline (§|H1). 

Figure 1131 redisplays the same full model set as in 
Figs. 1101 and llll but now color coded by the total burned 
mass. Models with higher Mbum (high Ek) generally oc- 
cupy the faster declining portion of the plot while models 
with lower Mburn have overly broad light curves. It is 
clear from the figure that a large part of the spread in 
the model WLR is due to the variations in Mburn among 
the models. Eliminating both the very high and very low 
Mburn points does result in a loose inverse correlation be- 



tween peak brightness and decHne rate (i.e., brighter is 
general broader), though there is still a large dispersion 
compared with observations. The best agreement with 
observations is achieved if common SN la burn 1.0 to 1.2 
Mq of their mass to siUcon and heavier elements. With- 
out further modification, models that burn 0.7, 0.8, and 
1.3 M0 seem to be rare events. 

We can further quantify the effect of Mburn on the 
light curves by examining a subset of models with fixed 
Mni = 0.5 M0 and Mpo = 0.1 Mq, but with Mime var- 
ied from 0.1 to 0.7 Mq. The total burned mass among 
these models thus varies from Mburn = 0.7-1.3 Mq cor- 
responding to a variation of Ek from 0.68 to 1.43 B. 
Figure^lshows that the models obey an anti-Phillips re- 
lation. Models with higher Mburn have shorter rise times 
and faster declines, and are also brighter at peak. Quan- 
titatively, increasing the amount of burned IME mass 
by 0.2 Mq (a Ek increase of 0.25 B) increases Mb by 
0.1 mag, AMi5{B) by 0.2 mag, and decreases the B-band 
rise time by 3 days. 

6.2.2. Effect of the Mass of stable iron and Mixing 

The second important model parameter leading to the 
dispersion in Figs. 1101 andllll is the mass of stable iron, 
Mpe. Since it is not a source of radioactive decay energy, 
and since its opacity is not so different from ^^Ni, the 
chief effect of Mpc is to influence the distribution of ^^Ni. 
For models with larger values of central Afpc, the ^^Ni 
center of mass is pushed farther out, thus decreasing the 
</)Ni parameter (this is not necessarily the case for ^^Fe 
and ^^Ni produced in the ^^Ni zones because of a finite 
metallicity and neutron excess). This can be expected 
to lead to faster light curves for three reasons discussed 
in §13 (1) based upon Equation the diffusion time to 
the ejecta surface should be shorter; (2) the occurrence 
of ^^Ni at lower densities leads to a lower percentage of 
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Fig. 15. — Effect of changing the mass of stable iron-group 
elements (M-pe) on the _B-band light curve for a subset of M-series 
models. Same as Fig. 1141 expect that here the mass of ^^Ni and 
the explosion energy are held constant (Mbum = 1-1 ^^0i in a-H 
cases), while the amount of stable iron is varied at the expense of 
IME. Here MNi= 0.5 Mq, Mp^= [0.0,0.1,0.2,0.3 Mq] and Mime + 
Mpe = 0.6 Mq. Models with larger Mpg have 56Ni distributed 
closer to the surface, and hence faster light curves. A change in 
Mpo of only 0.1 Mq gives the entire width of the observed band. 
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Fig. 16. — Effect of the degree of "mixing" on the B-band light 
curve. Same as Fig. 1141 but for four models all with the same 
ejected composition, Mni = 0.5 Mq, Mp^ = 0.1 Mq, Mime = 
0.5 Mq, but each with different mixing prescriptions. The more 
heavily mixed models have ^*Ni distributed closer to the surface 
and hence faster light curves. 



gamma-ray trapping in the post-maximum epochs; and 
(3) the increase in iron group elements at higher veloc- 
ity layers of ejecta leads to the stronger development of 
Fe II/Co II features in the post- maximum spectra, con- 
tributing to a faster color evolution (and hence B-band 
decline rate). 

To demonstrate the important effect of A/po in Fig. El 
models were selected with fixed M-^i— 0.5 M0 and fixed 
Mburn= 1-1 Mq, but with Mpc varied from 0.0 to 0.3 Mq. 
The models also obey an anti-Phillips relation - for given 
Mni, models with larger Mpo have faster light curves. 
Quantitatively, increasing Mpc by 0.1 Mq increases Mb 
by 0.05 mag, AMi5{B) by 0.1 mag, and decreases the 
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Fig. 17. — Same as Fig. 1101 but including only those models 
which obey the observational constraints of rapid expansion and 
"reasonable" nucleosynthesis and spectra, M^um = 1.1 it 0.1 Mq, 
A^IME > 0.1 Mq, Afpo = 0.1 - 0.3 Mq. 



i3-band rise time by 2 days. 

These effects of Mpo on the light curves are not re- 
lated to the presence of stable iron group elements per 
se, but to the effect Mpo has on the distribution of ^^Ni. 
Essentially the same effect can be demonstrated more di- 
rectly by varying the degree of mixing in the model. In 
Fig. 1161 we show four models each with the same com- 
positional production (Mm = 0.5 Mq, Mpe = 0.1 Mq, 
.^^IME = 0.5 Mq) but each with different degrees of mix- 
ing. The more heavily mixed models have a greater 
proportion of ^^Ni in the outer layers of ejecta, and 
hence faster light curves, for the same three reasons given 
above. 

7. TOWARDS A WORKING MODEL 

Having identified the two principal physical parame- 
ters that cause dispersion in the model WLR - the ex- 
plosion energy and the distribution of ^^Ni in the ejecta, 
selected subsets of the models are now examined that 
are in better accordance with the observations. What 
are the common properties of those "viable" models that 
fall within the observed WLR in Figs. EH and and 
why is the observed scatter so small? 

7.1. Constraints from Nucleosynthesis and Nuclear 
Physics 

First, one might consider what is "reasonable" from 
other quarters. Nature does not just select, with equal 
frequency, any random combination of nuclear products 
that unbind the star, but must select values consistent 
with known nuclear physics and, on the average, with the 
requirements of stellar nucleosynthesis and SN la spec- 
troscopy (e.g., the presence of IME in the spectrum). 
Does the application of these observational constraints 
serve to increase the agreement of the models with the 
WLR? 

It is not reasonable, for example, that a SN la make 
no stable iron, but only ^^Ni. Ignition in models near 
the Chandrasekhar mass can only be achieved at densi- 
ties in excess of about 2.5 x 10^ g cm""^ and any burning 
near that density makes ^'^^^^Fe and 58,60]sjj^ ^^Ni. 
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Fig. 18. — Same as Fig. 1121 but including only those models 
which obey the observational constraints of rapid expansion and 
"reasonable" nucleosynthesis and spectra, Mturn = 1.1 ±0.1 Mq, 



IME 



> 0.1 Mq, M| 



iPc = u.i — U.J iviQ. Once again, the solid 
line is the observed relation offPhiTl ips et al. (1999), with a cali- 
bration Mb = -19.3 for AMi5{B) = 1.1. The systemmatic offset 
of ~ 0.06 magnitudes in color between the observations and mod- 
els suggests that assumed thermalization parameter etjj = 1, likely 
overestimates the actual redistribution probability. 
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Fig. 19. — Left: WLR for a more finely gridded set of the mildly 
mixed (M-series) of models which obey the constraints Mtum = 
1.0 - 1.1 Mq, Mpe = 0.1 - 0.3 Mq, and A/jme > 0.1 Mq. The 
^^Ni mass varies fro m A f^jj = 0.4 — 0.8 Mq, with the same color 
coding as in 

Figure Right: Same as left, but including only 
those models that have total iron group production in the range 

A/pe + A/n; = 0.7 - 0.9 Mq. 



Further, if the star has any appreciable metaUicity, the 
excess neutrons will m ostly end up in these same stable 
iron-group isotopes. As iTimmes "etaD (|2003^ discuss, all 
initial CNO will end up at the end of helium burning 
in the isotope ^^Ne, creating a mass fraction approxi- 
mately 0.02 (Z/Zq). Subsequent burning and conser- 
vation of neutrons turns this chiefly into ^^Fe and ^®Ni 
with a combined mass fraction « 0.05 Z/Z©. For a range 
of ^^Ni masses up to 1 Mq, and metallicities Z = to 
3 Zq, this implies a stable iron mass of to 0.15 Mq. 
An additional minimum of 0.1 Mq of ^''Fe is expected 
from electron capture. This is roughly the amount of 
burning required to reduce the white dwarf density be- 
low the point where electron capture is important; (e.g., 
iNomoto. Thielemann, fc Yokoil [T984D . Thus we expect 
MpeS. 0.1 Mq always. 

On the other hand, nucleosynthesis in the Milky Way 
Gaiaxj requires that the sum of ^''Fe and ^*Ni be approx- 
imately 10% by mass that of ^^Fe. Given that massive 
stars also make some iron, this might possibly be raised 
to ~15%, but the iron in Type II supernovae is made in 
a region of high neutron excess as well. This suggests 
that Galactic SN la do not, on the average, make more 
than 0.2 Mq of stable iron per event. Of course, one 
does not know the isotopic composition of iron in distant 
galaxies. It is reasonable, however, that whatever phys- 
ical constraints operate to limit the amount of electron 
capture in local SN la also function in similar explosions 
far away. In summary, it seems that the stable iron mass 
is restricted by nuclear physics and nucleosynthesis to 
typically 0.1 to 0.3 Mq even for metallicities as high as 
three times solar. This does not mean there cannot occa- 
sionally be supernovae with very different characteristics. 
Nucleosynthesis constraints only operate on the average. 

The total mass burned is also limited by observational 
constraints on the typical expansion speed. The super- 
novae do not come apart with only a small excess of total 



energy over the binding energy, or spectral lines would be 
too narrow and ionization stages too neutral. In terms 
of light curves, it is already known from Fig. 1131 that ap- 
proximately 1.1 Mq needs to burn. This means that the 
typical model burns most of its mass. That is natural in 
all detonation models, and also true of strong deflagra- 
tions. 

Finally, the presence of strong silicon, sulfur and cal- 
cium lines in the maximum-light spectrum of SNe la im- 
plies the production of at least 0.1 Mq and probably 
0.2 Mq of IME for normal events. For these elements, 
there is no nucleosynthetic upper bound and the spec- 
troscopic limits are not presently highly constraining. 

Figure El a-nd Figure El show the resulting plots of 
Mbvs AAfi5(i?) and Mb vs B-V when these constraints 
are applied. In particular, we choose those models with 

Mburn = 1.1 ± 0.1 Mq, 0.1 > Mpe < 0.3 MQaud 

Mime > 0.1 Mq. Unlike Fig. EH the more restricted 
models do show a WLR, albeit a noisy one. The scat- 
ter in the color-brightness plot is also reduced. Nature 
apparently realizes this limited set of models more fre- 
quently. The WLR exists, not because a single parame- 
ter, the ^®Ni abundance determines both the brightness 
and decline rate of SN la, but because of other physics 
that constrains the production and distribution of stable 
iron, ^^m, and IME. 

The correlation between color and peak magnitude 
continues to be tight in Fig. El though the models are 
offset to the red about 0.06 magnitudes. This again prob- 
ably reflects a thermalization efficiency of less than 100% 

(§EU. 

7.2. A Constant Mass of Iron Group Elements 

The scatter of model points in Fig. El is still much 
greater than observed. To reduce it further, more strin- 
gent restrictions must be applied. In order to fur- 
ther populate the allowed band, addition M-series mod- 
els were constructed obeying the following constraints: 
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Fig. 20. — Effect of varying A/ni while keeping the total iron 
group production (Mt-ji + M-pe) fixed. The models shown have 
AfNi = [0.5,0.6,0.7,0.8] Mq, Mime = 0.3 Mq, and Mp^ + Mni = 
0.8 M0. Such a variation might be attributed to a variable amount 
of electron capture and metallicity in an otherwise standard super- 
nova. However, Model M080003, which has no neutronized iron, is 
not realistic. The total range of metallicity and electron-capture 
effects is probably bounded by the red and green points, i.e., about 
0.3 magnitudes. The range of metallicity effects alone is less. 



Fig. 22. — Effect of varying the abundance of stable iron group 
throughout the ^^Ni zone, as might be attributable to variations in 
progenitor metallicity. The figure shows six models derived from 
Model M070103, in which the ratio of stable iron to ^®Ni is varied 
throughout the ^®Ni zone from to 25%. All models have, in 
addition, 0.1 Mq of stable iron at the center of the ejecta. The 
total mass of stable iron thus varies from 0.1 to 0.275 Mq. The 
models follow the observed WLR. 




Days Since Explosion 

Fig. 21. — Same as Fig. [221 but computed with STELLA. A 
WLR of the correct sign is also present, though the agreement 
with observations for these particular models is not as good. For 
models with a small quantity of initial iron, especially M080003 
and M070103, the light curve at early times is fainter in STELLA 
because of a deficiency of lines included for Ni and Co (the Fe line 
list is more nearly complete). See text. 



Mburn = 1-0 - 1.1 Mq, Mpo = 0.1 - 0.3 Mq and 
Mime > 0.1 M©. A smaller range of '"^^Ni masses were 
also explored Mm = 0.4-0.8 Mq. Figure [T^fleft panel) 
shows the width-luminosity relation for this set of mod- 
els. The scatter is reduced compared to the full plot, but 
is still larger than the observed. 

A second cut of these models is then made assuming 
that the total mass of iron group elements, as well as 
the total burned mass, must also lie within a restricted 
range. Figure [T^ fright panel) shows only those models 
in which Afpc + A/ni — 0.7 — 0.9 Mq, which are seen to 
be the models that fall within the observed WLR. 



The necessity of having iron group elements in the re- 
gion 0.7-0.9 Mq follows from the spec troscopic and color 
evolution effect discussed in IKasen fc Wooslcv ( 2006) . 
The i?-band decline rate is largely determined by the 
rate at which Fe 11/ Go II line blanketing develops in the 
post-maximum spectra. In order for this line blanket- 
ing to develop, there must be a high abundance of iron 
group elements out to layers v ~ 8000 km s~^ which cor- 
responds to mass coordinate m = 0.8 Mq. This can also 
be achieved by mixing. 

It is not unreasonable that the explosion produce a 
nearly constant mass of iron group elements. The iso- 
topic composition may vary because of ignition density 
and metallicity, but all matter that burns with a temper- 
ature over ~ 5 X 10^ K will be iron-peak isotopes of some 
variety. Apparently, not only do common SN la burn a 
nearly constant fraction of their total mass, but a nearly 
constant fraction of that achieves nuclear statistical equi- 
librium. Delayed detonation at a nearly constant density 
(after a nearly constant amount of mass has been burned) 
might be one way, but not the only way of achieving that 
requirement. 

7.3. Variable Electron Capture and Metallicity 

The results of the last section suggest that SNe la mod- 
els in which both the total iron group production and the 
total mass burned are constants of the explosion are in 
better accord with the observed WLR. In such a scenario, 
the amount of Mni is varied at the expense of stable iron 
group elements such as ^''Fe and ^*Ni. There are two 
ways in which this may come about. 

First, stable iron group elements are produced at the 
center of the ejecta from electron capture (embodied in 
our Mpc parameter). If the central density of the white 
dwarf is higher, electron capture will be enhanced and 
the ratio of Afpc to Mni will increase. Such variations in 
ignition density might reflect variations in accretion rate 
in binaries with different separation, masses, metallicity, 
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etc. A second mechanism that will vary the ratio Mni 
to M-pc is the progenitor metallicity, which determines 
the relative abundance of ^^Ni to ^^Fe when burni ng to 
nuclear statistical equilibrium ijTimmes et al.ll2'003j) . 

These two effects are difficult to separate in the model, 
though ther e are observational con straints on the effect of 
metallicity ijGallaeher et al.ll2005(l . In the models, elec- 
tron capture produces a more centrally concentrated dis- 
tribution of stable iron, while increased metallicity makes 
the iron in the same place as the ^^Ni. But the effect on 
the light curve of increased electron capture plus exten- 
sive mixing may be indistingu ishable from that of in- 
creased metallicity (though see lHcifiich et al.lfl998L mix- 
ing all the way to the surface has observational diagnos- 
tics. We have in mind mixing that does not extend to 
the outer layers.). 

These two effects are demonstrated in Figs. 1201 1^ 
and 03 Figure EOl shows M -series models in which AfNi 
is varied from 0.5 to 0.8 M©, but Mpe + Mm = 0.8 M©, 
and Mime = 0.1 Mq are both held constant. These 
models, in which the stable iron remains concentrated at 
the ejecta center (as expected from electron capture) are 
in reasonably good agreement with the observed WLR. 

Even better agreement is found if the stable iron is 
mixed throughout the ^^Ni zone, as expected from vari- 
ations in metallicity. Figure |221 shows a set of models 
derived from M070103 in which a constant ratio of ^''Fe 
to ^^Ni is varied from to 25% throughout the ^^Ni zone. 
All models have, in addition, 0.1 Mq of stable iron lo- 
cated at the ejecta center. The total stable iron mass 
in these models thus varies from 0.1 to 0.275 Mq. The 
models fall very nicely along the observed WLR. 

Thus variations in either (or, more likely, both) the 
progenitor metallicity and the degree of electron capture 
can be a significant source of the observed luminosity 
variations in SNe la. However, it appears that this effect 
has a limited range. Varying the metallicity from to 3 
times solar only leads to variations of 15% in Mni- Larger 
values of metallicity may be unreasonable, and a floor, 
~ 0.1 Mq, on the lowest ^^Fe abundance is set by elec- 
tron capture (e.g., Nomoto. Thielemann. & Yokoi 1984). 
Less than 0.3 magnitudes of the observed peak magni- 
tudes, and perhaps 0.4 magnit udes in decline ra te can be 
explained this way (see also H ofiich et al1ll998fl . This is 
roughly one-third of the observed spread in peak magni- 
tude for common SN l a, and is con sistent with the obser- 
vational limits of Gallagher et al . ( 2005) that metallicity 
is not the major cause of the width-luminosity relation, 
or of the preponder ance of bright SN la in late galaxies 
ijHamuv et al.lll99'fi^ . This is especially true, given that 
part of the spread in Figs. 1201 and 1^ is due to electron 
capture during the explosion. 

Figurel^presents the same series of models as Fig. 1201 
but computed with Stella. Reasonably good agree- 
ment persists for models that keep both the total iron 
group production (MNi+Afpc) and IME production fixed, 
though the agreement is less good for brighter, broader 
supernovae because these models have no or little Fe in 
the beginning. Although Stella has ~ 10^ lines of Fe 
in the opacity, the line list is probably too poor for Ni 
and Co. Lower opacity means lower emission according 
to Kirchhoff 's law (which is applicable here when the line 
opacity is treated as an absorptive one). This makes the 
brightest models under luminous. 
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Fig. 23. — Effect of varying A/jMi at the expense of IME in moder- 
ately mixed models (Fig.0. The parameters are similar to Fig. 1141 
except that the mass of stable iron is held fixed while the mass of 
IME varies: Mni = [0.4,0.5,0.6,0.7,0.8,0.9] Mq, Mp^ = 0.1 Mq, 
and Mime + AfNi + Afpe = 1.1 Mq. The S-band peak magnitude 
declines even faster with decreasing AMis{B) than the observa- 
tions show. This is because moderately mixed models with large 
masses of IME and low ^^Ni have their ^^Ni more centrally con- 
centrated and thus dec line more slowly. Better results are obtained 
with more mixing ('Fig |25l . 



We note that our results in this section differ from 
the conclusions of IMazzali &: PodsiadlowskH lj2006|) , who 
also studied the effect of varying the ratio of stable iron 
group species to ^''Ni in SN la models. We find that 
varying this ratio leads to models in general a ccord with 
the WLR, whereas LMazzaU & Podsi adlowskil 12006) find 
that the variation creates dispersion from it. They there- 
fore identify the ratio as a possible "second parameter" in 
the WLR. The different conclusions likely stem from the 
different approa ches to the radiative transfer problem. 
IMazzali fc Pods iadlowski ( 2006) use monochromatic LC 
calculations and a simplified form of the mean opacity 
which depends only on the iron group abundance and 
time. These calculations therefore do not directly cap- 
ture the critical dependence of the model spectral/color 
evolution on the ejecta temperat ure and ionization state. 
IMazzali fc Podsiadlowskil l|200fiD do attempt to quantify 
the color effects using static synthetic spectra at select 
epochs. These spectral calculations, however, may be 
limited by the adoption of an extended inner boundary 
surface which emits blackbody radiation. In this case, 
the continuum formation is treated only approximately, 
and the radius of the inner boundary surface becomes an 
important free parameter of the calculation. 

7.4. Constant Mni + Mime 

Might the good agreement found in the last section 
persist if ^^Ni is varied at the expense, not of Mpc, but 
of IME? The models considered here have fixed AIp^— 
0.1 Mq and Mni varied from 0.4 Mq to 0.9 Mg at the ex- 
pe nse of IME. The procedu re is similar to that employed 
bv lPinto fc EastmanI l)2001|) . but not identical because of 
issues of mixing and, to a lesser extent, composition and 
explosion energy. The total iron group production here 
(Mpo + A^Ni) varied from 0.5 to 1.0 Mq while the burned 
mass was held fixed at Mbum= 1-1 Mq. 
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Fig. 24. — Percentage of trapping of gamma-rays from radioac- 
tive decay as a function of time for tlie models of Fig. 1231 in wliich 
M^i was varied. Models with lower Mni have ^^Ni concentrated 
closer to the ejecta center, and hence more efficient gamma-ray 
trapping in the post-maximum epochs. 



Fig. 25. — Effect of varying A/ni at the expense of IME, but in 
explosions in which the inner 0.8 Mq has been more thoroughly 
mixed (the "M C" Model in FiglTl. Unlike the stratified M-series 
models in Fig. I2ijl the MC-models better reproduce the observed 
WLR. 



Figure [221 shows that the slope of the WLR for this 
model set is steeper than observed. The slight positive 
correlation between Af^ and AAfi5(i?) is due to the de- 
pendence of th e spectroscopic /co l or evo lution on Mni, 
as discussed in lKasen fc Wooslevl l|200fiD and §|3 above. 
However, three other effects act counter to this in the 
moderately mixed models: 

First, in the low Af^i models, the distribution of ^^Ni 
is concentrated closer to the ejecta center. This leads to 
an increased effectiveness in the trapping of gamma-rays 
from radioactive decay, as shown in Fig. 1241 Near maxi- 
mum light {tcxp ~ 20 days) the differences in gamma-ray 
trapping among the models are minor, but by day -t-15 
{tcxp ~ 35 days) trapping in the A^Ni^ 0.3 Mq model is 
nearly 40% greater than that in the AfNi= 0.9 Mq model. 
This extra trapping slows the decline rate in the lower 
AfNi models. 

Second, again because ^^Ni is more highly concentrated 
towards the center, the optical depth (and hence diffusion 
time) to the ejecta surface is significantly larger in the 
low AfNi rnodels (see Eq. . This further tends to slow 
the decline rate of these models. 

Finally, as mentioned several times already, the B- 
band decline rate is determined largely by color evolu- 
tion controlled by the development of Fe II and Co II 
lines. In the low Mm models, the iron-group elements 
extend to only very low velocities. For example in the 
AfNi = 0.3 Mq model, the edge of the ^^Ni zone is at 
5000 km s~^ compared to 8500 km s~^ for the AfNi= 
0.7 Mq model. This inhibits the development of strong 
Fe II/Co II line blanketing in the low Mm models, and 
hence slows the decline rate of these models. 

The very steep WLR seen in Fig. reflects the fact 
that these three effects act counter to and nearly cancel 
the spectroscopic, color-evolution effect. Thus, in this 
particular set of calculations, we find that simply varying 
the ^^Ni mass at the expense of IME in a well stratified 
medium does not give a WLR relation with the correct 
slope. Either this represents a failure of the assumptions 
underlying the radiative transfer calculations, or (as we 




O.a 1.0 1.2 1.4 1.6 1.8 2.0 
4M,5{B) 




1.0 1.2 1.4 1.5 1.8 2.0 
4M,5{B) 



Fig. 26. — Left: WLR relation for the entire set of mixed core 
(MC-series) models. These models ar e ch aracterized by heavy mix- 
ing below m = 0.8 Mq. As in Figure lT^ they obey the constraints 
Mburn = 1.0-1.1 Mq, Mpg = 0.1-0.3 Mq, and Mime > 0.1 Mq. 
The ^^Ni mass varies fro m M j^; = 0.4 — 0.8 Mq, with the same 
color coding as in Figure ESI Right: Same as left, but including 
only those models that have total iron group production between 
J\/fc + Mni = 0.7 - 0.9 Mq. 



discuss presently) at least a moderate degree of mixing 
of the zones appears necessary. 

7.5. Centrally Concentrated Mixing 

The models in § 17.41 failed to reproduce the observed 
WLR because the distribution of iron group elements 
was more centrally concentrated for the low Mm models. 
This can be rectified if the the inner regions of ejecta 
are heavily mixed, thus equalizing the distribution (but 
not quantity) of ^^Ni among all models. Such mixing 
could, for example, reflect the maximum extent of the 
deflagration region in a delayed detonation scenario. The 
outer part, having been burned by a blast wave, would 
be less heterogeneous. If the density were low enough. 



Type la Supernova Light Curves 



17 



that part would be predominantly IME. ^ 

Observationally, the absence of significant high veloc- 
ity iron in the spectrum suggests that vigorous mixing 
was restricted to regions sO.8 M© (Branch ct al. 2005) . 
Transitions to detonation, if they happen, are believed to 
occur when the burning enters the "distributed regime" 
at de nsities ^ 1 — 3 x 10'' g cm~^ (jNiemever fc WooslevI 
Il997|) . Depending upon the detailed flame model and 
especially how the burning is ignited, such densities are 
reached wh en the fl ame has m oved through '^l Mq (e.g., 
lNomoto~^h ieleni ann. fc Yokoi 1984). This motivates 
the studies of models where the mixing has not been uni- 
form, but is concentrated in approximately the inner 0.8 
M0 of the star where it was assumed to be severe. The 
MC-series of models was prepared for a restricted range 
of burned mass, 1.0 and 1.1 Mq; ^^Ni masses, 0.4 to 
0.9 M0 and stable iron masses, - 0.3 M0 and a flner 
sampling of the intervals was also used. These models 
had the same flnal elemental yields and explosion ener- 
gies as the corresponding models in the standard "M" 
series, but the composition was distributed in a different 
way. These are called the "MC" (for "mixed core" ) mod- 
els. Figure^shows an example that should be compared 
with the corresponding "M" model. 

In Fig. |55l we show the equivalent of the models in 
Fig. 1231 but for the MC series. The mixing of the inner 
layers leads to an increased agreement with the observed 
WLR. Meanwhile, Fig. [211 (left panel) shows the full set 
of MC-series models. Compared to the equivalent set of 
M-series models (Fig.^1 left panel) the mixing serves to 
signiflcantly reduce the scatter and improve the agree- 
ment with the observed WLR. A more restrictive set of 
the MC models in which Fe + = 0.7 to 0.9 Mq fits 
the WLR even better (Fig. I^Hl right panel) . 

This flnal plot shows that it is possible to deflne a well 
defined, physically motivated cut of the model templates 
that does replicate the observed WLR. The chief req- 
uisites are: 1) a nearly constant explosion energy (and 
burned mass) around 1.0 - 1.2 Mq; 2) ^^Ni varied at the 
expense of either stable iron or IME, but sufficiently well 
mixed in the inner 0.8 Mq that its distribution in velocity 
space does not vary greatly, even when the composition 
is changed; and 3) a nearly constant iron plus ^^Ni mass 
around 0.8 ± 0.1 Mq. The resulting WLR is stiU quite 
broad though, filling the observed band with no excess 
left for observational errors and systematics, and there is 
no a priori reason why the explosions must always be so 
constrained. 

8. CONCLUSIONS 

The multi-band photometry (and, in some cases, the 
spectra) of a large set of Type la supernova models have 
been studied using two approaches to the radiation trans- 
port problem. These models were constructed in a sim- 
ple fashion (§ I2.1|l that allowed complete control over 
the major parameters affecting the outcome - the explo- 
sion energy and the masses of ^^Ni, stable iron, and IME 
produced. The models were exploded by depositing, uni- 

^ Detonation waves at intermediate densities (p ~ few X 10^ g 
cm""^) give explosive temperatures in the range 3.5 — 4.5 X 10^ K 
wliere burning produces a mixture of ^^Ni and IME. Lower densi- 
ties give only IME plus neon, sodium, and magnesium. Eventually, 
the density is too low (and the heat capacity of the radiation field 
too high) for any burning. 



formly, an amount of energy corresponding to the change 
in composition in a standard white dwarf model. All 
sensitivity to the uncertain physics of the actual burning 
was thus absorbed into the three model parameters, plus 
some prescription for mixing. Comparison of the light 
curves and spectra of sample models with those of more 
complicated models and to observations of SN la gave 
good agreement. 

Using these models, which consisted essentially of ev- 
ery explosion one could produce starting from a 1.38 
Mq carbon-oxygen white dwarf, the necessary conditions 
were determined for reproducing the observed relation 
between decline rate and peak magnitude (both in the 
i3-band). The physical parameters leading to intrinsic 
dispersion in that relation were also identified and quan- 
tified. The relation between peak magnitude and color 
at peak magnitude was studied as well, and found to be 
more robust to parameter variation. 

The set of all explosions with positive kinetic energy 
did not give a WLR. Instead a large range of decline 
rates, quite inconsistent with observations, was found for 
SN la with the same mass of ^^Ni. If SN la light curves 
are a single parameter family, that parameter cannot be 
as simple as just the mass of ^^Ni. Other possible cuts 
of the model set were thus considered based upon total 
kinetic energy (or mass burned), and various restrictions 
on the final composition and mixing. 

In order to satisfy the WLR, the ejecta structures of 
SNe la must obey certain constraints. First, all SNe la 
must have a common total mass burned (and hence ki- 
netic energy). In the present calculations, good agree- 
ment with observations was found if that common burned 
mass was 1.1 ± 0.1 Mq and the kinetic energy at infinity 
was 1.2±0.2 X 10^^ erg. Second, the radial distribution of 
iron group elements (including ^*^Ni) in the ejecta must 
be fairly uniform among all SNe la. We found the best 
agreement with the observed WLR among those models 
where the distribution of iron group elements extended 
from near the center of the ejecta to '^0.8 Mq. For ex- 
ample, one can consider the subset of models in which 
the total mass of iron-group elements is nearly constant 
at ~ 0.8 ±0.1 Mq, with variations in Mm arising from 
difference in the ratio of ^®Ni to stable iron group species 
produced in the explosion. Differences in this ratio are 
indeed predicted to arise from metallicity and electron- 
capture effects, although it is not clear that the entire 
range of SN la luminosities can be so explained given the 
constraints provided by the underlying explosion physics 
and Galactic nucleosynthetic measurements. Various ar- 
guments lead to the conclusion that the mass of stable 
iron (mostly 54,56pg ^^^j 58,60]\jj-j jg generally restricted 
to the range 0.1 to 0.3 Mq with about half of that range 
due to metallicity effects. Alternatively, one could vary 
the amount of ^^Ni at the expense of IME. In this case, 
good agreement with the observed WLR is found as long 
as the inner layers of ejecta (< 0.8 Mq) are fairly well 
mixed. In the present calculations, well-stratified ejecta 
structures do not reproduce the observed slope of the 
WLR, as the lower Mni models have ^^Ni concentrated 
relatively closer to the ejecta center, which tends to slow 
their light curve decline rate. To populate the entire ob- 
served span of the WLR, especially at faint luminosities, 
the mass of these IME is quite large in some instances, 
arguing for a delayed detonation. 



18 



Woosley et al. 



The present calculations adopted an LTE approxima- 
tion for both the level populations and the line source 
functions (i.e., the thermalization parameter, eth = !)• 
Although a common approach in many previous trans- 
fer studies, this (and other simplifying approximations 
adopted) can be expected to lead to errors in the model 
observables on the order of 0.1 to 0.3 mag. The rela- 
tive differences between models, however, may be known 
more accurately, and so the general trends obeyed by 
sets of models and the level of dispersion in such trends 
may be considered more robust. For example, our choice 
eth = 1 overestimates the true redistribution probabili- 
ties and leads to maximum light colors that are too red 
by ~ 0.06 mag relative to the observations in the color- 
peak-magnitude plot of, e.g.. Fig. ^] Nevertheless, the 
clustering of the points in this plot is quite tight and 
the slope is correct. The possible refined calibration of 
the radiative transfer calculations and its impact on the 
study of the WLR will be explored further in a subse- 
quent paper. 

Given the large dispersion found, even in selected sub- 
sets of our models (e.g.. Fig. [T7jl . the narrow spread in 
the observed relation between Mb and AAfi5(_B) is sur- 
prising. No single physical parameter yields it, unless 
several other parameters are highly constrained. Part of 
that spread can be reduced by a judicious choice of mix- 
ing and a tight constraint on the mass burned (Fis.H^. 
Even so, it is difficult to avoid the conclusion that a major 
fraction of the observed scatter of the WLR, ftO.l mag, 
reflects intrinsic physical diversity, and not observational 
effects. This conclusion has important implications as 
one plans for future studies using SN la as calibrated 
standard candles to ever higher precision. Given the dis- 
persion expected from the models themselves, a larger 
sample of supernovae will have to be studied to get a 
strong signal to noise ratio for cosmological effects. 

A primary benefit of the approach adopted in this 



paper is the ability to construct an expansive grid of 
SNe la explosion models without restricting ourselves to 
any specific theoretical explosion paradigm. Such a large 
and general database of model light curves and spectra 
nicely compliments the large sample of SNe la data sets 
currently being acquired by ongoing observational pro- 
grams. A direct comparison of the models to individual 
observations would be useful in interpreting the phys- 
ical properties of any given SN la event. In addition, 
the model grid should be helpful in developing refined 
techniques for calibrating the luminosities of SNe la so 
as to limit and reduce intrinsic sources of dispersion or 
evolution. For example, our studies suggest that an im- 
proved strategy might be spectroscopic template fitting 
in which as much data as possible from each individual 
SN la is compared to a template of model spectra and 
colors. The models employed could be the ones computed 
here or some derivative of that set. 
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